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a Abstract. We propose a class of random scale-free spatial networks with nested community struc- 

tures and analyze Reed-Frost epidemics with class related independent transmissions. We show that the 
epidemic threshold may be trivial or not depending on the relation among community sizes, distribution 
^— H of the number of communities and transmission rates. 

> 

C^ Keywords and phrases: Epidemics, SIR, Reed-Frost, percolation, long-range, directed, scale-free, modular, 

^^ nested, communities, hierarchical, threshold, spatial. 



f^ 2010 Mathematical Subject Classification: 



X 



Primary 60K35; 92D30. 
Secondary 93A13; 82B43; 92C60. 



1. Introduction 



We consider a spatial random graph which at the same time is scale-free and has a nested community 
structure, and study Reed- Frost SIR epidemic ([E], [25]) on it. We find that with a natural transmission 
mechanism, in which transmissions occur independently with rates related to community sizes, the critical 
threshold is trivial or not depending on the relation between community sizes, distribution of number of 
communities to which each individual belongs and rate of the decay of the transmission probability as the 
community size increases. Scale- free networks ([S], [2], [T3]) have been widely studied in the context of 
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epidemics (see [55] and [TU]) suggesting at first that this might lead to triviality of the critical threshold 
( |24) ■ |17) . [TS], [31] )■ On the other hand, most scale-free networks lack a spatial dimensionality, which is 
quite relevant to make the models more realistic (see e.g. [H]): one of the few prosed networks possessing 
both features has been suggested by Yukich 33 . Yukich's network is, however, missing network modularity, 
i.e. the gathering of individuals in communities with faster transmission rates (see [7], [8], [3], [4] and [5]), 
a feature which has gained recent interest due to its relevance in infectious transmission. The formation of 
communities can be described by several mechanisms, such as random intersection in which extra vertices 
randomly connect to the vertices of the graph and links are then generated between vertices connected to a 
common extra vertex (see [11] for a description of random intersection and a review of other mechanisms) . 
However, most real community structures are nested (see, e.g., |30j . |32j and [12j ) unlike the networks 
generated by random intersection and similar mechanisms. 

The class of random networks discussed here have spatial features, are scale-free and possess a nested 
community structure. The networks are based on a connectivity graph, which, for simplicity, is here taken 
to be Tfi endowed with a hierarchical structure of partitions into larger and larger communities. To generate 
the network each vertex v ^Tfi \s assigned a random integer value X^, where the X„'s are i.i.d. random 
variables. Each vertex v identifies an individual, which belongs to all communities up to level X^ in the 
hierarchical structure. The basic random connectivity graph is obtained by adding to the nearest neighbor 
edges of Z'' all the edges between pairs of vertices belonging to at least one common community. For a wide 
class of distributions of the X„'s the connectivity graph is scale-free. 

We then consider Reed-Frost SIR epidemics on the connectivity network, in which infected individuals 
at time t contact each neighbor independently with some transmission probability, and if the neighbor 
is susceptible it becomes infected. To complete the model, it is natural to consider basic transmission 
probabilities for nearest neighbor vertices, and then an additional probability, decreasing with the size of 
the community, of independent transmission for any community shared by two individuals. In this way, the 
transmission probabilities do not depend only from the connectivity graph, but directly from the shared 
classes, and give rise to a very realistic mechanism. The set of individuals ultimately affected by the Reed- 
Frost SIR epidemic is the set of vertices belonging to a percolation graph with connection probabilities given 
by the transmission probabilities ( •21] ) ; for natural choices of the probabilities of infection through shared 
communities the phase diagram of the percolation graph exhibits a transition from nontrivial to trivial 
percolation threshold. 
In summary, the model depends on five parameters: 



(i, indicating space dimension; 

z, determining the growth factor z'^ of community sizes; 

a > 1, determining the distribution of the number of communities to which an individual belongs; 

p, indicating the transmission probabilities to neighbors; 

p, modulating the decrease in transmission probabilities for large communities. 



Several random networks can be generated along the indicated lines. In particular, the construction must 
specify the form of each partition and the interconnections between partitions. To illustrate the mathematical 
properties of the networks, we discuss in Section 2 a very simple and schematic structure, in which at each 
level fc the space is partitioned into hypercubes of linear size z*^, which are then packed into hypercubes 
of linear size z^^^ and so on. To keep things simple one can think to z = 2. For simplicity, wc also limit 



ourselves to just one single parameter a to generate the connectivity graph, although this is excessively 
simplified, as the inclusion in small communities is likely to follow a different pattern from that of inclusion 
in large communities. In Section 2 we give a detailed description of the construction of the connectivity 
network. 

In Section 3 we show that the degree distribution Dy of any vertex v in the connectivity network satisfies 



1 d'^/^+i 



^d 



2a ' d — log^ a 



lim P{Dy > h) h-^+^ e 

/i— >-oo 

where ojd is the volume of the d-dimensional unit ball and 7 — 1 = , °f - " , so that the network is scale-free 
for all a e (1, z^); in particular, for z^ < a < z^ the network exhibits the typical value of 7 e (2, 3). 



In Section 4 we complete the description of the Reed- Frost epidemic and begin the description of the phase 
diagram in the a — p variables; such description is completed in Section 5 by dominating the probability of 
transmission in certain sets by those in a long range percolation model extending a recent result in [22] ; it is 
remarkable that although we use edge variables to bound a model based on site variables the result is still 
sharp and we identify the exact phase diagram. 

(1) For a > z'^ the network has short range behavior, and the hierarchical communities structure is 
irrelevant for the existence of critical threshold: there is a critical epidemic threshold Pc for all p. 

(2) For 1 < a < z"^ the behavior depends on p: if p < -^ there is still a nontrivial epidemic critical 
threshold Pc G (0, 1) while Pc = for p > -%. This means that percolation, and thus an infinite 
outbreak, occurs at all values of p in the parameter range we just identified. In the scale-free region, 
determined by a G (l,z''), Pc is thus trivial or not depending on the transmission rate in large 
communities. It is trivial if the transmission rates are constant {p = 1) or with a not too fast rate; 
on the other hand it is not trivial for p below a critical curve in the phase diagram. 

(3) On the line a = 1 each vertex belongs to all communities: the model is similar to long-range 
percolation (see j27j ) and is studied in ^23 : Pc ^ or pc > for the same parameter range as in 
long-range percolation. 

In a sense the proposed model interpolates between short (a > z ) and long {a — 1) range percolation. A 
summary of the phase diagram is in figure 1. 



2. The connectivity graph 



Consider a random graph Ga.z = (Z'^, ^a,^) with Z*^ as set of vertices and a random set of edges Eq,^^ to be 
specified. In the first place, Ef C Ea^z, where Ef is the set of nearest neighbor edges of Z'^. Then, consider 
i.i.d. random variables Xy, v G Z'^, with a nonnegative integer distribution pct,v such that Pa,v{Xv > fc) = 
a"*^, k — 0, 1, ... , where a > 1 is a parameter. We let pa — Yivez-' Ma the joint product distribution of 
the Xy^s on the Borel cr-algebra in X = N^ . By this choice there is only one parameter determining the 
distribution of the number of communities to which individuals belong; the average number of communities to 
which an individual belongs, a measure called group membership (see [H] and [H]), is J2k '^'^ = i'^^ 1)^^- 
This is a realistic number especially for a G [5/4, 2). 



Next, let z > 2 be an integer and for each k partition Z'' into blocks 

Bz,k{i) = B^^kih,-- .,id) = {w = (^'i,-- ■,Vd) eZ"^ : z^i^ < Vj < (ij + l)z'' - 1, for all j ^l,...,d}. 

Blocks represent a system of nested communities. Note that vertices separated by coordinate hyperplanes 
lie always in different communities; the community structure is thus confined to orthants, and vertices in 
different orthants are connected only through nearest neighbor connections: this is not an unrealistic feature, 
however, as it might represent very rigid borders or seas. 

Given pa and the Bz^kiiYs, the random connectivity graph Ga,z is completed by including into the edge 
set Ea^z, next to the nearest neighbor edges, also all pairs {u,v} such that 3k G N,i £ Z"^ with X„,X„ > 
k and u, u € Bz^k{i)- In other words, given a, z, /i^ and Bz^kiiYs, the random graph Ga,z is defined by a 
map (j)z: X ^ H = {0, f }^^ where E'' = {{u,v} ■.u,ve 1'^}, with 

4>z{x){u,v} = I{3feeN, iel.'^\x^,x^>k and u,veB^_k{i)} ^ ^^{Itt— u|=l} 

by Pa,z = IJ'ai4'7^)- Later we are interested not only in the connectivity graph but also in the set of 
communities joining each pair of vertices: this leads to further specify the map 0^, as done in section 4 
below, but we first study the connectivity properties of G^.z- 

For V el.'^ and 1] e H, let Dy = Dy{ri) = \{u e T^ : {u, v] e Ea,z}\ be the degree of v. 
Lemma 2.1. For all v £l/ and a, z such that a < z'^ 

(1) lim P„,,p, >/i)-/i^-i>-^, 

h—too Lot 

where 

d - log^ a 

Proof. Given w e Z'' and /i e N consider the block Bz^i(h) = Bz^i(h){i) such that v G Bz^i(h) with l{h) ~ 
I -r^ log, h+l\. Then 

Ld— log, a oz J 



E^o 



By the CLT, lirnii IJ-aiY.ueB.^n^-, hx^>i{h)} >h)> 1/2. Hence, 

/iT-i Po,,p. >h) > /i-'-i ^lc.{X, > l{h)) P„,, f J2 h{v.u}eE^_,} >h\X,> l{h) 



n 



To get the corresponding upper bound on the degree distribution we compare the connectivity graph 
to Yukich's network, which has vertex heights based on uniform distributions and connections related to 



distance. As first step, we compare the connectivity graph to a network based only on distances but retaining 
the distribution of the Z„'s for the vertex heights: for (5 > consider G'^ z s ~ C^"^' ^L z s) such that 

(2) {u,v) e E'^.^g <:=> 3fc s.t. Xu.Xy > k and d{u,v) < Sz''; 
more precisely, let (j)'^ g : X ^ H such that 

(3) 4'z,si^){u,v} = I{3fceN I X„,X„>fe and d{u,v)<Sz''}{^) 

and let P'^ ^ g — fXai{4'z s)^^)- Note that for S = Vd and for every increasing A C H 

(4) P^^"(^)<IPUy5(^)- 

In fact, taking k = min(X„,X„), if d{u,v) > \fd z^ then {u,w} ^ £^q,z- Therefore, if {u,i'} G £^q,z then 
z*^ > d{u,v)/y/d, so that {u,v} G E' ^, i.e. (?!)z(a;) < </>' /;j(2:)- This implies that HA is increasing, a; G X 

and (j)z{x) G A then also 0^ /^(a^) € ^, i-e. ((^z)~^(^) C (^'^ ^j (A). 
Note also that 

(5) <l>'z,six){u,v} ^l^^x^^z^^> <'(Y) y{x). 

3. Comparison with Yukich network 

Let {Uv}v,zz,d be i.i.d. uniform [0, 1] random variables with distribution Pjj on [0, 1]^ and consider Yukich 
network Gs.s — ^'^ ■, Es,&) defined for s, (5 > by 

(6) {u,v} G ^,,5 ^ d{u,v) < 6mm{U-',U-') 
As before one can take W = [0, 1]^ , define (j)'^ g : W -^ H such that 

and let P'^, ^ = Pu{{4''s s)~^)- We need to slightly reformulate Theorem 1.1 in Yukich ([33]) to incorporate 
the constant S. 

Proposition 3.1. For all d,6 and s G (3,00) 

8 hm t^ KADAv) >t)=[ -T-^ 

t->oo ' \ sd — 1 

for all u G Z"^, where LOd denotes the volume of the unit hall in W^ . 



Proof. Yukich proves the same result for 5—1. The conclusion is achieved by taking the origin v = Q, 
conditioning ow Uq — t and using translation invariance. The basis of Yukich proof is Lemma 2.1 in |33) . 
which states that 

E(i^,(0)|[/o = t)k f \x\-i dx = duJd r t^-^-i dt = -^p^ r-t^-^-i) 

J\x\<T-= Jo sd~l 



When a generic 6 is considered we get 



E(D,(0)|C/o = r) w / \x\-^St^ dx 

J\x\<5t- = 



dujdSi I t''-'--i dt = :i4:rr:i^ r 



sd-1 



,d I 1 1 'y— 1 



m — z^ 



The rest of the proof in ^331 is still valid with the constant /3 = f^zf replaced by ^"f/^^^ U 

Then we can deduce the following upper bound for the power law distribution of the network Ga,z- 
Theorem 3.1. For all z and a £ (1, z'') 

(9) lim PaAD{v) > h) h^'^ < (- 
where 7 - 1 = d°tg% 

Proof. First recall that 

(10) P.,aP>/i)<P:^„_y5p>M 

since {D > h} is increasing. We want to compare G'^ z 5 to the Yukich's network G'^ g. For 

l^a [z^" > m) = /i„ (X^ > fc) = a"*^ = a" '°s. ™ ^ ^^ '°s. « 

On the other hand Pu{U~'^ >m)~ Pu{Uv < ?n^^/*) — m^^/" so that taking log^ a — 1/s we have 

PuiUy" > m) = ^aiz^" > m) for m ^ z^ 

and 

Pc/ ([/„"' > m) = a-'°g^™ > ar-i°s='"T = ^,(X„ > [log, ml) = ^^^(X^ > log, m) = ^^^(z^" > m) 

for all other m e M. Therefore, C/i, and z-^^ can be coupled by the following joint distribution. Let P[/„ be 
the distribution of U^ and v^ be a probability on the cr-algebra in [0, 1] x N such that for A C [0, 1] and 
fc e N it holds 

uM: k) - Pu. (a n (^-('=+l)/^ z-^-/^]) . 

We have 

. vM.^) ^ PuAA); 

. z.,([0, l],k) = FyJlz-e^+DA-.z-W/^]) = a-'^ - a-C^+D = ^^^(X. = k); 

• Vy{{u,k) : u^"" > z''} = i^y{{u,k) : u < z^''^''} = 1. 

The product distributions fia and Pjj can be coupled by the product distribution v — YiveZ'' ^"j under which 

f/^" > z^" for all V G Z"^ with probability one. 

If w e W^ and x G N^ are such that w^'^ > z^" for all v then 






Thus, if A C iy is increasing, x e (0'^ ^) ^{A) and w^^ > z^^ then w G ((/)(, ^) ^{A). Hence, for A increasing 

Pi^AA) = M« (('/''.,.)-' (^)) 
= v(w,{ct>'^A-\A)) 
= ^{{4>\io^^.)-^^)-\A),{cl^'^A-\A)) 

Since A — {D > h} is increasing 

P. AD >h)< Pi^^^/D >h)< P(;,^^„)_.^^p > h). 

If we take s = (log^ a)^^ and a G (1, z'^) then s e (1/d, oo) and the resuh follows from Proposition 3.1 with 
S = Vd. D 



From lemma 



2.1 



and theorem 



3.1 



for large h it holds P^^aiD — h) k, h ^ where 7 — 1 = j^°f^ "^ . 



-logj 

Thus the hierarchical model is scale free for each a G [1, z'^). Typically, in the scale free region —3 < —7 < —2, 
which is then equivalent to z^s < a < z^ . 

We end this section by commenting on the relation between the scale free region and the average number 
of communities to which an individual belongs. As we have seen, there is a realistic average number of 
communities for a £ [5/4, 2], which has no intersection with the typical scale-free region even for d = 2 and 
z = 2. It is, however, quite simple to realign the parameter ranges by introducing some more parameters more 
realistically describing small group membership. This is reminiscent of long range percolation in dimension 
1, in which the probability of nearest neighbor connection can, by itself, determine phase transition for a 
critical value of the main parameter ([!]). We do not pursue this direction here. 



4. Epidemics 



We consider a Reed-Frost dynamics to describe the spread of an infection on the connectivity network (see, 
for instance, [11], section 3, for a detailed description). In such dynamics, at discrete times each infected 
individual contacts each one of its neighbors with some probability, and if the neighbor is susceptible it 
becomes infected; in the meantime the infected recovers. Differently from usual, we assume, however, 
that the probability of infectious contact depends on the communities shared by the two neighbors: in 
particular, we assume that there is a probability of independent transmission for each community shared by 
two individuals, and we are interested in the set of individual eventually affected by the epidemics started 
from one single vertex, the origin for instance. Such set can be identified with the cluster Vg containing 
the origin in an edge percolation process on Ga,z described by the following probability measure: for each 
value a; e M'* of the Xy's, consider a (conditional) Bernoulli probability distribution Fx,z,p,p on {0, 1}^"='^ 
such that 

00 

(11) ^x,z,p,p{V{u,v} = 1) = 1 ^ _[ _[(! ^ W ^{{it,"}! SieZ'': x„,x„>fc and u,v<£Bk,^ii)}) ■ 

fc=0 



Our main interest here is in studying for which values of the parameters there is a finite or an infinite set 
of infected individuals, or, equivalently, a finite or infinite cluster, i.e. we are interested in the probability 
]P£c,z,p,p(|Vo I — oo). The joint probability distribution which describes percolation and epidemics is defined 
on the Borel u-algebra in H by Va,z,p,p = Jx ^x,z,p,p ^laidx). 

For a given x let G^z be the realization of the connectivity graph with value x of the X^ variables. Since 

Gx,z contains all the nearest neighbor edges and they are open with probability at least p, if p > ttc , 
the critical point for d-dimensional bond percolation, then percolation occurs regardless of the value of the 
other parameters and of the realization x. Notice that ttc < 1 by Peierls argument, and, more precisely, 
TTc = 1/2 ([16 ) and ttc ^ l/2d ([20 ). Moreover, for any fixed x and p, the probability in ( 11 ) is increasing 
in p, and the random variables 'q{u,v} ai'c independent. Thus, it follows by a standard FKG inequality (see, 
e.g., [in]) that for any p > p' and any increasing event A C {0, 1}^"^.== we have Vx,z.p,p{A) > ]Px,2,p,p'(^)- 
Since ^q oo = {^ ■ ^o I = °°} i^ increasing it follows that there exists a critical Pc{x, p) < 1 for the onset of 
an infinite percolation cluster. 

It could happen that Pc(x,p) — 0. If a = 1 then we are assuming Xy = oo for all v, and the percolation 
model is quite close to long range percolation ([1?]) in which the critical threshold tTc has a transition at 
some value of a parameter which corresponds to p: for small values of p we have tTc > and for large p it is 
instead tTc = 0. After showing that, in fact, the critical threshold Pc{x, p) = is almost surely constant in x, 
we see that a similar transition occurs in the hierarchical model for all values of a G [1, z"^]. To this purpose 
we introduce a more detailed description of the model: consider E — N^ x {0, 1}"^ ^^ and parameters a, p 
and p. Then take a Bernoulli probability distribution Pq,^,^ on the Borel cr-algebra ^ in E such that 

• Pa,p,p{(^v >k)= a^'^ for all v € Z''; 

• Pa,p,p{<^{u,v}.k = 1) = PP'' for all {u,v} e E'' and fc G N 

One then retrieves the probability Pa,z,p,p by considering the map ipz '■ "S -^ H such that 

and observing that Pa,z,p,p = Pa,p,p{'4'z^)- Notice that while Pa,p,p is Bernoulli, the distribution Pa,z,p,p on 
H is not independent since, for instance, if m,m' G -Bz,fc(0) \ i?z_fc_i(0) then Pa,z,p,p{V{o,u'} — l|?7{o,u} = 1) = 



a 



-fc -I ^,-2A; 



^ a ^^ = Pa,z,p,p{il{o.u'} = 1)- Pa,z,p,p IS actually one-dependent, xxxx 



Lemma 4.1. pc(x,p) is almost everywhere constant in x 



Proof. Under Pa,p.p the variables (t„'s and cr(/„ „} ^j's are collectively independent. Consider the a-algebra 
An generated by the variables with index in {{v, {v, u}, k) : v & {[—n, nj'^nZ'^y; {u, v} G E'', u, w G ([— n, n]''n 
Z'^Y; k > n}; then Aoo — <^nAn is trivial under Pa,p,p- 

Since the event A^ = {r]\3v G Z'' : \V} ■^ | = oo in 77} is such that ■!/)J^(Aoo) G An for all n G N, then 
il;~^{Aoo) G Aoo and Pa, p,p{ip7^ i^oo)) — 0, 1. Thus, Aoo has probability zero or one for Pa,z,p,p-8:a.. rj £ H. 
Hence, Px,z,p,p{Aoo) =0,1 for /^a-a.a. x £ X. Since Pc{x,p) exists for all a; G X, it is pa almost surely 
constant. D 



We can define Pc{ct, p) = ini{p : Pa,p,p{'4'z ^(^oo)) — !}• We already know that pc{ct,p) < 1. We see now 
that pc{a,p) — when the transmission probabihties for large communities do not decrease fast enough. 

Lemma 4.2. For a £ [1, z'^) and p > -^ we have Pc = Pc(a, p) — 0. 



Proof. The joint probability Pq,p,p suggests several dynamic constructions of the epidemics together with the 
reference graph; one is the following. Starting from the origin consider the sequence of boxes Bz^k = ^z,fc(0), 
fc = 1, . . . and sequentially generate the following variables: 

(0) ^o; 

(la) a-u.v e Bzs; 
(Ifc) (^{o,v}s,v e Bz.i; 

(ka) (7v,V e Bz^k \BzM-i; 

ih) <^{v,u},j,u,v G Bz^k~i \ Bz,k^2,j = 1,. .. ,fc - 1; 

(fcc) a-{v^u},k,u € Bz,k-i,v e Bz^k', 

(Last) (yiu.v},o for s-U nearest neighbor pairs {u, v}. 

Note that at every step only new a variables are generated, that the last step can be performed at any time, 
possibly subdivided in several steps, and that the procedure generates all relevant a variables in the positive 
orthant: in fact, if u',u' G Bz^k \ ^z,fc-i then criv^ujj is generated at step {{{k + !)(,) for j = 1, . . . , fc and 
(jc) for ah J > fc + 1; if, instead, v' e Bz^ \ Bz.k-i and v' G Bz,k+r \ Bz,k+r-i, r > 1, then cr{„,„}j for 
j = 1, . . . , fc + r — 1 is not generated but it is also not relevant in the process and for j > fc + r is generated 
at step (jc). 

Following this construction we can show that for a G [l,z''), p > -% and any p > there is an infinite 
cluster. We generate a sequence ik,k G N, of vertices in Bz.k \ -B^.fe-i or empty sets with the following 
procedure, in which the definition of ik depends on 3 events which may occur depending on the status of 
ik-i'- 



• if ctq ^ 1 then io = 0, else ig 



if ik-i G Bz,k-i \ Bz^k-2 and 3v G Bz.k \ Bz.k-i such that cr^, > fc + 1 and (T{i^_^.v},k = 1 then ik 

equals one of such vertices v (the first in some fixed order) ; 

if ife_i G Bz,k~i\Bz,k-2 and 3v G Bz,k\Bz,k-i : tr^ > fc + 1 but for aU such u's a[i^_^^y}^k = then 

ik equals one of vertices v with the first two properties (the first in some fixed order); 

if ife-i G Bz.k~i \ Bz.k-2 and for all v G Bz.k \ Bz^k-i we have a^ < k + 1 then ik = 0; 

if ik-i = and 3v G Bz.k \ Bz^k-i '■ cr^ > k + 1 then ik equals one of such vertices v (the first in 

some fixed order); 

if ik-i ~ and for all v G Bz,k \ Bz.k-i we have cr„ < fc + 1 then ik — 0; 



Given the vertices z^'s we can define the events: 



• 
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• Ak = {3v e B^^k \ Bz^k-i ■■(Jv>k + 1, cr{ifc_i,t,},fc = 1} 

• C'fc = {3u G Bz^k \ Bz^k-i : CTt, > fc + 1 but either ik-i = or for all such u's <J{i^_^.v},k — 0} 

• Ek — {for all w G -Bz,fe \ B^^k-i it holds (Tt, < fc + 1} 

where clearly Ak is not defined if ik^i ~ 0- Notice that all the events Ak,Ck and Ek are defined in terms 
of the variables at steps (ka) and (fee) of the construction outlined above. This implies that such events 
are defined in terms of variables which, once ik-i is given, are independent from those involved in defining 
Ai, Ci and £'i for i = 1, . . . , A; — 1. Moreover, for each k the three events form a partition of the probability 
space. Therefore, the sequence Zk ~ ak{ck,ek respectively ) if Ak{Ck,Ek respectively ) occurs, is a (non- 
homogeneous) Markov chain, whose transition matrix can be estimated in terms of the a variables. In 
fact, 

(12) P{Zk = ak\Zk-i = ak-i) = 1 - (l - ^) >!-« ^^' 

(13) PiZk ^ Ck\Zk-i ^ ek-i) = l-(l-^^^^j >l-e ^^' 

ppfc(^dfc_^d(fc-i)) 
and all other conditional probabilities are smaller than e 0*=+^ if Zk-i = ak-i or Zk-i = Ck-i and 

smaller than e =''+1 if Zk-i — Ck-i- 

We have 

P{Zk = ek)= 22 P{Zk^ek\Zk-i^z) P{Zk-i^z) < e ^^^i 

z=afc-l,Cfe_i,efc_i 

and 

P{Zk = Ck)= Y. P{Zk=Ck\Zk-i=z) P{Zk-i=z)<2e ^^^ 

SO that if p > a/z'^ 

00 00 

^p(Zfe = cfe) < 00, y^p{Zk = Ck) < 00. 
fe=i fc=i 

By the first Borel-Cantelli Lemma Ek and Ck occur only a finite number of times, so that with probability 
one the sequence terminates with one Ck and then A/j for h > k. In such case the vertex ik is connected 
to an infinite cluster containing all vertices i^ for h > k. Since there are countably many vertices there 
must be one k and one vertex v G Bz.k \ Bz^k-i which is starting vertex of an infinite cluster using edges in 
communities at level at least k with probability ci > 0. Such vertex can be connected to the origin using 
nearest neighbor edges, which are independent from the previous construction as they were involved only in 
the last step of the dynamic joint generation of graph and epidemic, with some probability C2 > 0. In the 
end, the probability of percolation from the origin is at least C1C2 > 0. D 



5. Domination by long-range percolation 

The description of the a — p phase diagram is completed by the following result. 
Theorem 5.1. For a> z'^ or a E [1, z''] and p < a/z'^ we have pc > 0. 
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This amounts to prove that, with the parameters a and p in the indicated region, there exists p > such 
that percolation does not occur for that value of p > 0. To show this, we actually bound the probability 
of existence of an infinite percolation cluster or infinite infected area in the nested model with that in a 
long-range percolation, for which it is easy to show that percolation does not occur for some values of the 
parameter by bounding it with a subcritical Galton- Watson process. 
A long-range percolation model is defined as a probability on the Borel cr-algebra in H such that Q^ s(j]{u v\ = 

^1 ~ {d(u.v)y 

Theorem 5.2. When s = \og^{a/p) and 13' = ^(i^)h^°s.d^ -^ ^^^^^ ^^^^ 



The main difficulty lies in the fact that in the nested hierarchical model the distribution on the edges is 
one dependent: we face this problem later on. Initially, we once again compare the percolation network to 
G^ 2 5 endowed with slightly larger infection probabilities than in the nested model. 
Let u,v e Z'', define 

kij(u,v) = llog^ — - — I, 

and consider a Bernoulli probability distribution P'a.z,p,p,s on the Borel a-algebra ^ in E' = M^ x {0, 1}'^" 
such that 



^'a,p,p,^ (K >k) = a-'' for all v e Z^; 

^'a!p!p,5 (aj„_,} = 1) = T^ pkiAu.v) fo^ all {u, v} e E^. 



Consider then the map ip'^ g : T.' ^ ¥,''■ such that 

i'^z,si'^')){u.v} = I{<,<>fci(M,«); ctJ^_^j=1)} 

Lemma 5.1. For all increasing events AC- H, f'a,p,piipjl{A)) < F'a,z,p,p,siiip'z s)^^i^)) 



Proof. Consider the cr-algebra Ax generated by the variables (T„,u € Z'', and let Vx^p.p = Pa,p.p{'\x) and 
P'x,p,p,(5 be the conditional probabilities of Pa,p,p and P'a,z,p,p,<5, respectively, given Ax- Notice that the 
conditional probabilities no longer depend on a and that F^^p^piip^^) ^^^ ^'x,z,p,p,s{{'^'z s)~^) are Bernoulli 
distributions on (the Borel cr-algebra of) H under which 

p..p,p(^."'('^{«M=i)) =1- n (i-w') 

where 

Ix{u,v) ~ {k\ 3? G Z ■.UjVE Bk^z{i) and Xu,x^ > k}, 
and 

P 



~^' .,z,p.:p.8 ((^L.)-'K„,.} - 1)) = T^ P'' 
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Note also that Ix{u,v) C {k-^ yd("> ^)i ^i JJii''^^'^) + 1; • • • i min(a;„, a;„)} so that 

1- n (i-p/'') ^ 1- n (1 



-PP') 



k>k^ /^(Uj-u) h>k>k-^ /^(ti,!;) 

since the series in the second Unc is alternating with decreasing coefficients. Therefore, 

P,,p,p(C'('^{n,.} = 1)) < P\,p,p,Vd iii^T\^{uM = 1)) 

and P'^ n p Vd ((V-'z)~^) dominates in the FKG sense Fx,p,p{ip'^'^)- 
Therefore, if A C H is increasing then 

p„,p,p(V;^(A)) = / p,,p,p(V'-^(A))/z„(dx)< / r,,,.p,p.Vd(«,y5)"'(^))A*o(^^)-P'a,p.p,Vd((^i)"'(^)) 

n 



To compare the percolation network P'^ ^ Jd^ii'' /ji) ^) with a long-range percolation network we are 
going to prove an analogue of Theorem 3.1 in [52]. In this direction there are two main problems. On one 
side, [22] applies to directed paths; on the other side, connectivities in [22] are described by convex functions 
k{Xy, Xu) and for values of Xu = Xu the connectivities are bounded by expected values Xy = E{k{Xy,Xu))- 
In that paper the reason why the connections become independent in different directions is that the Xy 's are 
constant. 

The directionality of the paths is easy to fix: paths under Pa z p p ^^^ '^ot directed, but can trivially 
be considered so by fixing an order along each path. Paths are instead ordered under P' ^ y/dii'^'z)~^) 
since the involved edge variables are defined according to an order. Theorem 3.1 of [Hj applies to hoppable 
collections of paths, such as the collection of all self-avoiding paths starting at the origin and reaching the 
boundary of some fixed set; since from each path one can extract a self-avoiding one. Theorem 3.1 applies 
to the occurrence of a connection from the origin to the boundary as well. 

As to the connectivity functions, the analogous in the present context would be k(Xy, X^) = i4'ziX))[y^u} 
which is not convex and cannot be easily related to any constant value. To proceed, we introduce families 
of i.i.d. random variables, one family for each v E Z"^, of the form X'/ ,,u G Z"^ \ {v}, and then bound 

P'^ z rt p \/d((^z)~"'^) by a network based on the X'/ .^s. Connections in different directions are independent 
and depend only on distances, thus the network based on X'/ ^n's is actually a long-range percolation model. 
This is possible if we take the probability that X'/ y^\ > k greater than or equal to the square root of the 
probability that X^ > k. This, in turn, implies that in the long-range model the presence of a vertex is 
equivalent, in distribution, to the fact that X^ > k for one of its end-points, say the smallest in some fixed 
order. While this implies that the probability that the infection travels a self-avoiding path is larger in the 



long-range model. Theorem 5.4 below shows the same inequality holds for the probability that at least one 



paths is travelled among those in a fixed suitable collection. 
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Consider (3 > and a Bernoulli probability distribution V" /3^z,p,p,s on the Borel cr-algebra A in S" 
N^'x^'\{('^*)''ezn X {0, ijS'' such that 



^~"p,z.p,pA<^'{^,v} >k)^ /3-'= for all {u, v)eZ''xZ''\ {{i, i),i e Z'^}; 
'~"p,..p,pA'^'(^,^,} = 1) = ^^4^ foi- all {"'^'} e E'^. 



Ph{<Tv>k) = 




Ph{(J{v,u) >k) = 


1 


^h(0-{«,n} = 1) = 


1- 



Consider then the map V'" 5 • 5^" — >■ -ff such that 

and let P^'^.^^^p^^ = P"/?..,p,p,5((V'"5)^'). We denote by P".",.,p,p,5((V'"5)"') the conditional probability 
given x" € X" = N^ ^^ \{(i,i),iei. } p^Q^g ^-j^g^^- jj^ passing from V'a,z.p,p,s to P"^,z,p,p,5 we have changed the 
network mechanism and kept the same transmission rates. 

We introduce an interpolation between P'q, z^p^p,^ and V" p^z,p,p,s- To this purpose we select an ordering of 
Z"* = {«!, W2, • ■ • } and, for /i = 0, 1, . . . , we consider the sequence of sets V^(0) = 0, . . . , V{h) = {vi, . . . , Vh}- 
For later purposes we take the order such that V{n'^) = B^ = [0,n — 1]'' n Z'^. Then we take a sequence of 
Bernoulh distributions Ph defined on the Borel cr-algebras Aih) of Y,{h) = N^^^('') x NV'(h)xz^{(^,^),^ez<^} ^ 
{0,1}^" by 

^ veZ'^\ vih) 

V e v{h),ueZ'^\v 

p 

Furthermore, define the map ''pz,h '■ '^{h) ^ H given by 

i'^Pz.5,h{(^)){u,v} = I{o-t(„,„)>fei,i(«,t;), crt(„,,)>/ci,5(«,t;), <T{„.„}=1} 

where i(w,i;) = u ii u e Z'^\V{h) and t{u,v) = {u,v) Hue V{h). We have Poiipj},) = Pa,^,p,p,^((V'^)"^)• 
Fix now a box i?„ = [0,n — 1]'' n Z'^ and consider the variables <j\b„, which are the u's restricted to 
Bn, i.e. to the index set {{v),{v,u),{v,u} : v,u e Bn}. For v,u e Bn and h < n'', {ip'z s hi'^")){v,u} and 
{'>Pz,5{cr)){v,u} depend only from ct"|b„ and a\B„, respectively. Therefore, P„d(V';^;^J = P" i3,z,p,p,s{itp")^}) 
by the definition of Pj^. 

Given a box Bn C Z'' and ?; G Bn, let i?i,_„ = {{w, u} : u e Bn D Z"^} and consider now a pair of (possibly 
empty) sets A, B C iJ„„, which in our case coincides with both E'^ and E* of [22], any |yl| -dimensional 
vector X — {xi, . . . ,x^a\) G (M+)I'^I and any |_B| -dimensional vector y — (j/i, . . . ,y\B\) G (M^)'"^'. For a fixed 
h, the values x and y are interpreted as realizations of X^ if u G V{h) or X(„ „) if w ^ ^(^)5 respectively. 

For A C £"„ „ we indicate by Za the event {77 : 77(1,, u} = for all {v, u} G A} C iJ that none of the edges 
of A is open, and for any probability P on H we define the zero functions Zy{P; n; A, B; x, y) — P{Za U Zb) 
as the probability that either none of the edges of A is open or none of the edges of B is open; for any 
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pair of probabilities P^"'' and p('') denote by z„(p(''),n) < z^{P^''\n) the fact that z^{P^°-^;n; A, B;x,y) < 
Zy{P^''^; n; A, B; x, y) for ah pairs of disjoint and possibly empty sets of endpoints A, _B C E^^n, all x e M'"*' 
and y G M'^L The extension of Theorem 3.1 in [22] that we are going to prove uses the following inequality. 

Theorem 5.3. If (3^ = a then for alln, /i e N such thatvh € Bn, z^^{Ph^i{'4i~\ /i-i); '^) — ^vh{Ph{^^\ /J'"-)- 



Proof. For fixed _B„ C Z'* and v ~ Vh & Bn, notice that the events ip^ I ^(^a) and ip^ ] u^^b) are measurable 
with respect to the variables a^, iT(„^„) and cr^^^^^^ which are indexed in the set Z„_„ — {v} U {{v,u},u e 
Bn \ {v}} U {{v, u), u& Bn\ {v}} . Then let A,B C £;„ „, disjoint, with |^| = r and \B\ = m, and x G M'^I 
and y G M'-^' be fixed; we identify each edge in A or B by its endpont different from v. We then let 

(14) AUB ^ {ui,U2,...,Um+r) 

indicate the vertices wihch are endpoints (different from v) of edges in AU B, ordered according to the 
distance of the endpoint from v, which is d{v, Ui) < d{v, Ui+i). We also indicate A — {vi, V2, ■ ■ ■ , Vr} and 
B — {wi, W2, • • ■ 7 Wm}- For simplicity of notation denote by d^ — d(y, Ui) the distance from v to Ui and by 
Oui , l^ui the following probabilities 

(15) a„. = /ia(^.>log.^)=a-'°^^(w) 



/3u. = /i,(x.,.. >log,^)=/3-'°^^(^) 



Thus a„. = (/3uJ^. Furthermore, let g^^ — ^^-^^ — ^ and Pi = Ph-i and P2 = Ph] we want to prove that 

(16) Pi(Z^UZb)>P2(Z^UZb). 

Let's proceed by induction on the cardinality of ^ and B. Note that if |A| = or \B\ — then Pi(Z^UZb) = 
V2{Za^Zb) = 1. 

(i) Suppose A = {u}, B = {w}. By symmetry we can assume that d^ < du] then a^ > Oiu and /3u, > /?„. 
We have 

= 1 - OiuQuqw 

¥2{ZaUZb) = i-P2(zinz|) 

= 1-P2(Z1)P2(Z|) 
= 1 - PuQu l^wQw 

Since /S^, > jS^ then /3„/3^, > /3^ = «„ and 

Pi{ZaUZb)>V2{Za^Zb) if|A| = |B| = l. 
In particular, equality holds if d^ — dyj. 



(ii) Now consider {ui, U2, ■ ■ ■ Um+r} — {^1,^2, • ■ • Vr} U {wi,W2, ■ ■ ■ Wm} — AU B such that d^ < du2 < 
■ • < o^Mm+r- Note that for any probability P 

?A u Zb) = HZa) + nzB) - F{Za n Zb) 
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As before, consider the probability of ZyiUZfi. With respect to Pi, if X^, < log^ -|J- then Z^i occurs. Instead, 

if fogj, -^ < Xjj < log^ ~^ then there exist i connections in the basic graph and Za occurs if at least one 
of them is open. Thus 

r—l j r 



j—1 2—1 i— 1 

m—l j m 



r+771— 1 



i^l 



i^l 



m+r 



J — 1 i— 1 i— 1 

With respect to P2, since edges are open independently of each other, we have 

r 

V2{Za) = X{{i~Pv^qv^ 

711 

1=1 

r m 

P2{ZAnZB) - Y[{l-Pv,qv,)Y[{l-l3^,qw,)- 

i=l i=l 



We proceed by induction on m + r: we show that if ( 16 1 holds for m + r —1 then it holds also for m + r. The 
vertex Um+r can be either in A or in _B and we assume with no loss of generality that Um+r — Vr (z A. Then 
we show that if Pi(Za' U Zb) >P2(Zyi'UZB) with \A'\ =r-l, \B\ = m then Pi(Za U Zs) >V2{ZaUZb) 
with A ^ A' U {v} and thus \A\ = r, \B\ = m. This is equivalent to show that 

(17) Pi{Za U Zb) - Pi{Za' U Zb) > P2{Za U Zb) - P2{Za' U Zb). 

By elementary calculation it turns out that 

r-l 



h{ZA) - Pi{Za') - -a^„qv„ Y[i^ - 9. J 



i=l 
r-l 



h{ZA n Zb) - Pi(Z^, n Zb) = -a„^g„^ H^l -'?".) IT (^ " '^-.) 

j=i 



4=1 



thus 



Pi(Z^ U Zb) - Pi{Za' U Zb) = -a.^q,,. n(l - Iv.) 1 - l[i^ - Iw,) 

i=i L j=i 

Similarly, with respect to P2 we have 

r-l 

P2{Za) - P2{Za') = -f3,,.qv,. Y[{1~ f3v,qvj 
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r~l 



i{ZA n Zb) - V2{Za' n Zb) = -/?.,g., Y[ii- /3.,g.J n(i - Z^-.?-.: 



j=i 



so that 



r-l 



j(Za U Zb) - ¥2{Za' U Zb) = -Pv^v^ \{i^^ ^-.9">) 1 " H^l " Z?-,?-,) 



i=i 



In order to prove inequality (17) we must show 

r— 1 r rti -1 r— 1 r ni 

Since a^^ — /3^^ , this is equivalent to show that 

r— 1 r m -ir — 1 rm 

/?«. 11(1 -'?-•) i-n(i-9-.) <n(i-'^-.9-.) 1-11(1-/^-.'?'' 
?'— 1 J— 1 j— 1 

Since /S^^ < 1, 1 — qwi < 1 ^ l^wiQwi, thus 



j=i 



l[a~qv,)<Ylil-Pv,qv,) 



i=l 



(18) 

Moreover, we see now that 

m -| r 771 

(19) Pv^ l-[](l-g,„J < l-n(l-/3-,9« 

'- i=i -' '- i=i 

proceeding by induction on m. If tti = 1 then 

Pv^ Qw — Pw Qw 

because Vr is the vertex at maximal distance from u, so that (3^^ < P^- Next we evaluate the increment 
between the (m — l)-th and ?Ti-th term. 

rn — 1 n ra—1 

3 = 1 



/3.. 



1-11(1 -9-.^ 



/3.„ 



J = l -' '- J = l 



7TT. — 1 



i=i 

m— 1 

Pw^qw^ n (1 ~ Pw,qw,) 

j=i 



thus inequality ( 19 I follows from inequality ( 18 ) and {3^^^ < /3. 



n 



Now we are able to follow Meester and Trapman's work [22] to bound from above the probability of 
large outbreak, i.e. the existence of an infinite open path, by the corresponding quantity in the long-range 
model. In order to prove the results below we need to recall some definitions; the detailed definitions are in 
|22) . An ordered set of edges in some £' C Z'' x Z^ of the form S, = (wo'yi,wiW2, • ■ • T^n-iVn) is a (directed) 
path from vq to u„. A path ^ = (t;oWi,fiW27 • ■ • ,''^n-ii'n, ■ • • ) with infinitely many different edges is an 
infinite path. Given a finite or infinite path ^ = (wofi, Wif2, ■ • • ,^'n-ii'n) we indicate the truncation after 
k edges as ^"(fc) = (vt)Vi,viV2, ■ ■ ■ , Vk-iVk) and the tail starting after k edges as (,*{k) = {vkVk+i, ■■■)', for 
two paths fi — (wo'i^ij''^i'i^2j ■ • ■ , 'i^n-i'i^n) ^ud ^2 = (I'n'^n+i, • • • ) wc dcuotc the conjunction by {S,i,£,2) = 
{vqVi, 'yiW2, ■ . . , Wn-iWnj VnVn+1, . ■ .)■ Ncxt, let S be a collection of paths; if _£*•"-' is the collection of the first 
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n edges of E according to some given enumeration of E then we indicate by S„ the set of finite paths of S 
all of which edges are in E''"^ together with all the infinite paths of S truncated at the first instance they 
leave £'("1 

Furthermore, given a configuration rj E H = {0, 1}^ we say that ^ is open in rj if for all edges {vk, Wfe+i} 
we have r]{vk,vk+i} = 1- And we indicate by C" the event that at least one path in S is open. We say that 
S is hoppable if 

• for any v £ Z"^ and any two paths ^ and (j> of 'E. going through v, where v is the end vertex of the 
i-th edge of ^ and the starting vertex of the j-th edge of (j), then (^*(i)i '/'*(j)) € S. 

• lim„ C"" = C" 

Theorem 5.4. For every hoppable collection of paths S in E*^ 

(20) P'o.,.,p,,Ai^'.,s)-'ic-)) < p'^.^^.,.,.M:^s)-\c~)) 

Proof. We mimic the proof of Theorem 3.1 of [52], dividing the argument into 3 steps. Since P'a,z,p,p,s and 
P"r- r are not defined on the same space, we use the interpolating distributions Ph, which are such that 

two consecutive ones differ only in the variables related to a single vertex. Fix a box i?„ — [— n, n]'^ D U^ . 

(i) The first step is to show that for all 71 and /i such that w/i G i?„, P/i_i(V'j5 /i_i(C"")) < Ph{i'Z\ hk^^'^y) ■ 



Since/32 = a, byTheorem[5Jzi,^(P,,_i,n) <z^^(Ph,n). Denote by S'(/i) ^ N^ \v {h) ^^{v {h)\vu)y.-L \{{^.^),^& }x 
{0,1}^'\^-" = N^'\(^('')u-h) X N(v(/^-l))xz^{(^,^),^ez''} ^ {Q, 1}E'\s..- by T.'^(h) its restriction to B„, and 
J^y^ and A^ „ the Borel cr-algebras generated by the variables in S'(/i) and S^(/i) respectively. For all h 



PhX^;XkiC^'')\^'^'^ih))dPh-ii^kw). 

where for a' E E^(/i), P{ |cr') is the conditional probability given A'f^ ^; the last equality holds since i\ 
coincides with Ph-i on A'/^^ ,^. Therefore, 

Phii^lUiC-)) - Pk-ii^lU-riC-)) 



Now one can follow the proof of Theorem 3.1 in [22]: if the event C"" occurs in a'^, c^x regardless of the 
variables in Zvf^^m then the integrand is 0. Otherwise, one can follow verbatim case 3. of the proof of Theorem 
3.1 in [21 to conclude that Ph-ii4'~ih-AC~"K'jh))) < -P^l^.J, J^""ks;,(/,))) for all /i = 0, ... ,n'' - 1 
and thus the unconditional inequality holds. 

(m) By iteration, 

P'c.,..p,pAi^'.,sr\c~-)) = Po{^7loic-)) < p„.(CL-(C"")) = ^~"a,..p,p,5((C5)-'(c^-"))- 
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(Hi) In the last step we consider a general hoppable collection of paths ,=,. By definition of hoppable 
collection of paths, since C"" is decreasing in n, it follows that 



and using the previous steps the proof is completed. 



n 



Proof, (of Theorem 5.2). For all hoppable collections of paths .^, C" is an increasing event in H; moreover, 
{|Vq I = 00} ~ C" when S is the collection of all infinite paths containing the origin. If s = \og^{a/ p) and 
/3' = _P_(£^)^iog.<i^hen 



P (P 






1-p 



1 — p Va. 



idiu,v)Y 



Ql3',s{fl{u.v} = 1) 



for a long-range percolation model Qp'^s- Combining Lemma 15. 11 and Theorem 5.4 we have 



< Vo.,.,p,pAi^'.,6)-\\vr\=^)) 



< 



< 



(d)i 






D 



Proof, (of Theorem 5.1) In order to establish for which values of the parameters a,p^p,z no percolation 
occurs, it's now sufficient to dominate the long-range percolation model QjS'^s by a subcritical Galton Watson 
tree. Recall that a GW tree is subcritical, i.e. the probability of extinction is one, if the expected value of 
the descendants of any vertex is less or equals to one. If i?„ denotes the number of neighbors of a vertex v 
we have 






1 



P P' d(u,t;)'°s.(f) 

for all p e [0, 1] if a > z'^ or for a e [1, z"^] and p < -% 



<2dp + J22dk'^-'T^A-J-^ 



1 



ken 



1 



P P 



fJ°sAf) 



< 00 



n 
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Figure 1. The phase space of the nested model in the a — p plane. 
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